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Abstract 

We show that the imaginary part of the embedding potential, a 
generalised logarithmic derivative, defined over the interface between 
an electrical lead and some conductor, has orthogonal eigenfunctions 
which define conduction channels into and out of the lead. In the 
case of an infinitely extended interface we establish the relationship 
between these eigenfunctions and the Bloch states evaluated over the 
interface. Using the new channel functions, a well-known result for the 
total transmission through the conductor system is simply derived. 

1 Introduction 

In scattering theory, and consequently in the theory of electrical conductance 
PQ, the concept of channels is fundamental - channels are the asymptotic 
states at a particular energy, identified by some quantum numbers, between 
which the electrons are scattered. In conductance studies, which is our main 
concern, the current enters a nanostructure (for example) through open chan- 
nels in the leads on one side, and is transmitted into the open channels on 
the other side |2j. The theoretical problem in finding the conductance then 
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becomes one of calculating the transmission probability between the chan- 
nels J3J. In this paper we show that channels can be usefully defined in 
terms of the embedding potentials which embed the nanostructure on to the 
leads |H E] . Moreover, this gives a simple derivation of the well-known result 
[SJ El Ej that the total transmission through the nanostructure between the 
left- and right-hand leads at energy E is given by 

T lr (E) = ATriGUE^mE^G^E^mE^E)]. (1) 

Here Gi r is the Green function connecting the left- and right-hand surfaces 
of the nanostructure, and Ej/ r is a self-energy, or embedding potential, which 
couples each surface to the left- or right-hand contact respectively. The trace 
is over the quantum numbers of the channels in each contact. It was Wort- 
mann et al. jHj who realised that SmS;/ r in is the same as Inglesfield's 
embedding potential [I], and derived (JTJ) in the embedding context. They 
have used this result with embedding to study the transmission of metallic 
interfaces jH], and recently field emission has been calculated using the same 
formalism [Oj. 

The embedding potential in has been widely used in surface, interface, 
and transport calculations 01101 \U\ In surface calculations, for example, 
the surface region is joined on to a semi-infinite substrate, which can be 
replaced by an embedding potential. The embedding potential, added on 
to the Hamiltonian for the surface region, ensures that the surface wave- 
functions match onto the substrate wave-functions. This means that the 
whole problem can be treated by solving the Schrodinger equation for one 
or two layers of atoms at the surface rather than the whole semi-infinite 
substrate. 

In this paper we shall show how the embedding potential over a surface 
can be used to define channels crossing the surface. We shall relate these 
channels for the case of an infinitely extended two-dimensional surface of a 
semi-infinite solid to the Bloch states in the solid. Our results enable us to 
derive some interesting formulae in scattering theory, and in particular allow 
us to derive ((H) rather easily. 

We use atomic units (a.u.), in which e 2 = h = m e = 1, and also eV for 
energy (27.2 eV = 1 a.u.). 
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2 The embedding potential and flux 



The embedding potential E, which is the basis of this paper, is a generalised 
logarithmic derivative for the substrate region which it replaces [I], satisfying 

J^(rs) = -2 J s d\ s J:(r s ,r' s )^r' s ). (2) 

Here the integral is over the boundary S between the surface and substrate, 
and in the partial derivative ns is the normal to S, taken into the substrate; 
ip(r) is the solution of the Schrodinger equation in the substrate, integrated 
from some boundary value ip{rs) on S with outgoing boundary conditions 
into the substrate. £ is the same as the mathematicians' Dirichlet to Neu- 
mann map |12j . 

The embedding potential can be found from Go, the substrate Green 
function satisfying a zero-derivative boundary condition on S [A j , 

£(r 5 ,r' 5 ) = G -W^), (3) 

where Gq 1 is the inverse of Go over 5*. The surface inverse can be avoided by 
writing £ in an alternative form, in terms of a double normal derivative of 
the Green function G satisfying the zero amplitude boundary condition [T3] , 

1 f)G 

In this paper we shall be considering the embedding potential as a surface 
operator in the real-space representation of the Hamiltonian, but it can also 
be found in tight-binding form jSl E] , where it is usually called a self-energy 
The current can be written in terms of SmE, which is ultimately why 
this term appears in the transmission expression To show this, we start 
from the expression for the current density 

j = I{V>*vv-W}. (5) 

We now consider an embedding surface S, the interface, say, between the 
nanostructure and one of the leads (figure 1). Substituting (J2J) into (jSJ), the 
flux across S is given by 

I = i [ d 2 r s I d 2 r' s [r(vsMr s ,v f s )^(v' s )-^(v s )E*(v s ,v' s )r(r' s )]. (6) 

J S J s 

Note that because of the convention for the sign of the normal derivative in 
(|2*|). this is the flux from the nanostructure into the lead. As the embedding 
potential has the symmetry property of Green functions, namely 

£(r 5 ,r' s ) = £(r' 5 ,r s ), (7) 
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Figure 1: Nanostructure joined on to metallic leads, represented by shaded 
areas. Region I, consisting of the nanostructure + right-hand lead is embed- 
ded over S onto the left-hand lead, region II. Subsequently it is convenient 
to label the left- and right-hand interfaces by Si and S r . 



equation © simplifies to [S] 

I = -2 f d 2 r s f d 2 ^f(r s )SmS(r Sl r' s )i|i(ry. (8) 
J S J s 

The fact that (JBJ) gives the flux into the lead is associated with our convention 
that £ corresponds to outgoing waves. This is because G , whose surface 
inverse gives the embedding potential (JHJ), is evaluated at energy E + ie, 
where e is a positive infinitesimal. 

At this stage we introduce the eigenfunctions of 3mS at each energy as 
a natural choice of channel functions, given by the equation 

f rf 2 r^mE(r Sl r' s )^(r' s ) = A^(r 5 ). (9) 
J s 

Although Q only gives the channel function ipi(rs) over S, this uniquely 
defines the corresponding outgoing wave-function throughout the substrate. 
Because 3mS(r5', r' s ) is real, and symmetric in r$ and r' 5 , the eigenvalues 
are all real. The eigenfunctions can be taken to be real and orthogonal. If 
the eigenfunctions are normalised to unity over S, 

d 2 r s ^s) 2 = 1, (10) 



s 
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we see by substituting into (jSJ) that the flux associated with ipi is — 2Aj. 
This must be positive or zero, because of the outgoing boundary conditions 
implicit in the embedding potential, so we conclude that SmS is negative 
semi-definite, with eigenvalues negative or zero. The eigenfunctions with non- 
zero eigenvalue correspond to flux-carrying, open channels; there is an infinite 
number of closed channel eigenfunctions satisfying Q with zero eigenvalue. 

SmS can be expanded in the usual way, as a sum over its eigenfunctions, 
and with the normalisation given by (fT0|) we have 

3mE(r s , r' s ) = £ A^(r s )^(r' s ). (11) 

i 

This reduces to a sum over the open channel functions, 

3m£(r s ,r' 5 )= £ X^s^'s), (12) 

open i 

a result which we shall find useful later on. 

In conductance, it is only the open channels which are important, those 
with non-zero eigenvalue. The same open channel functions can be used for 
flux out of the lead into the nanostructure; this corresponds to using the £ 
appropriate to incoming waves, which just involves a change in the sign of 
SmE. These channels have the advantage of being defined via (JHJ) over the 
interface between the lead and the nanostructure, precisely where they are 
needed. They do in fact correspond to the usual channels in at least one 
case, the case of waveguides. With a waveguide lead, the separability of the 
Schrodinger equation for motion along the waveguide and across the guide 
leads to different waveguide modes, and an embedding potential which is 
diagonal in a mode representation [TT]. As a result, each embedding potential 
channel function ipi(rs) corresponds to a waveguide mode. 

3 Bloch states and channel functions 

In the case of an infinitely extended interface, it is interesting to study the 
relationship between our channel functions and the substrate Bloch states 
evaluated over this interface. With an extended interface, the "lead" on 
one side consists of a semi-infinite metal in which all wave-functions can be 
labelled by the Bloch wave- vector parallel to the interface, K. The embedding 
potential can be expanded as 

E(r 5 , r' s ) = S K,nmXK,n(rs)XK,m(rs) (13) 

K n.m 
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using a set of orthonormal surface functions XK,n( r s)> i n general complex, 
that satisfy Bloch's theorem with wave-vector K. Using (JJJ) it follows that 

3mS(r s ,r' s ) = ^S K (r s ,r' s ) (14) 

K 

with S K expanded as 

S K (r 5 , r' s ) = SK,nmXK,n(r S )XK, m ( r 's) ( 15 ) 



n,m 

where 

y -I 

2z 



1 V~l* 



(16) 



Sk is Hermitian and so has real eigenvalues, and its eigenfunctions are the 
channel functions since the flux across the surface S of an outgoing state 
with Bloch wave- vector K 

^( r s) = Z^ C K,nXK,n(r S ), (17) 

n 

is given by 

I = ~ 2 c K,n^K,nraCK,m- (18) 

At each K and energy E there are Bloch states travelling towards and 
away from the interface, as well as solutions decaying exponentially into the 
bulk ^SJ. These Bloch states and evanescent states can be considered as 
open and closed channels respectively, and channels defined in this way have 
been used in transport studies across interfaces [HI HE!- In line with our 
convention, let us take the open Bloch channels to be the states propagating 
away from the interface, into the bulk; the states travelling in the opposite 
direction simply involve the perpendicular component of the wave-vector 
changing sign. We shall denote these Bloch states and the exponentially 
decaying states by 0K,i( r )> with % labelling the states at fixed K and E. 

These Bloch states, evaluated over the interface, also diagonalize Sk 
(hence QmE), but we shall see that they are not the same as our channel 
functions ipK,i- As before, the flux associated with K ,i is given by 

h = -2 J s d 2 r s | s d 2 r^(r 5 )SK(r s ,r' s )0K,(r' s ). (19) 

There is also a well-known result 0, which follows from Green's theorem 
and the Bloch property of the wave-functions, that for i ^ j 



I d 2 r s 
Js 



W'*)-g->)-^(r5)-^(r*) 



0, (20) 
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hence 

/ d 2 r s [ d 2 r^(r 5 )S K (r 5 , r' s )0 Kli (r' 5 ) = 0, i ? j. (21) 
J S J s 

Unfortunately, equations p9j) and (|2*T|) do not imply that the 0k/s are eigen- 
functions of Sk or SmS. Although these functions are orthogonal over the 
three-dimensional unit cell, it is known that the 0k,i's are not orthogonal 
when integrated over the surface, 

/ dV 5 0^(r 5 )0 KlJ (rs) ^ 0, i± j. (22) 
J s 

However, the 0k,i's are linearly independent - this makes it possible to use 
them in surface and interface matching procedures. In order for the function 
space spanned by the open 0K,i( r s) and ipK.,i{ r s) to be the same, there must 
be the same number of Bloch functions as our channel functions carrying 
finite flux. 

The transformation matrix, which relates the flux-carrying Bloch states 
with the open channel functions, has the very useful property of being uni- 
tary. Let us normalise the open 4>K,i( r s) an d ipK,i( r s) to carry unit flux per 
two-dimensional unit cell, away from the interface, so that 

J d 2 r s J d 2 r' s <j>*^<j>^ = ~ 5 f, J d 2 r s J dV^E^ = ~y • (23) 

We now expand the flux-carrying Bloch functions 4>K,i, 4>kj in terms of the 
open and closed channel functions (the closed channel functions can be nor- 
malised to unity using (fTUJl ). 

0K,i(r,s) = Yl aim^K, m (r s ), 0K,i(r<?) = E ajn^pK,n( r s)- 

open + closed m open + closed n 

(24) 

Then from (|2*Hj) we obtain 

Jd'rsJd'r's^E^j^-^ 

£ a *m a jn J d 2 r s J d 2 r^^ m E K ^K,n- (25) 

open + closed m,n 

But the integral on the right of (f2*5*j) vanishes when m or n is a closed channel, 
and when m, n are both open the integral is given by the second equation in 
(J2HJ). Substituting, this gives 

Z) a *im a im = fy, (26) 

open m 
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Figure 2: Au(OOl) interface at K = 0: (a) band-structure in the TX direc- 
tion; (b) eigenvalues of £k=o on an embedding plane midway between atomic 
planes. The calculation uses a finite imaginary energy of 0.1 meV; (c) eigen- 
values of Ek=o on an embedding plane not midway between atomic planes. 
The numbers label the symmetry of the corresponding eigenfunctions. 

showing that the matrix a; m in the space of open functions is unitary. 

As an example, we consider the case of embedding onto semi-infinite 
Au(001), taking K = 0, zero wave- vector parallel to the interface. The band- 
structure as a function of k±, the perpendicular component of the wave- 
vector, is shown in figure 2a, corresponding to the TX direction in the bulk 
Brillouin zone. For comparison, the figure also shows the eigenvalues of Sk 
for K = as a function of energy, in figure 2b where the embedding surface 
S is midway between atomic planes, and in figure 2c where the embedding 
potential has been transferred to a plane that does not intersect muffin-tin 
regions of the Au potential (THHnj. We see that the number of bands is indeed 
the same as the number of non-zero eigenvalues, except at isolated energies 
where there are 5-function peaks. The 5-function peaks will be discussed 
later, but other features in the eigenvalue spectrum can be understood in 
terms of a one- dimensional nearly-free electron model. In this model the 
embedding potential near a band edge Eq behaves like [T%] 

3m£ oc -yJ\E-E Q \ (27) 
if the amplitude of the wave-function at Eq is finite on the interface. On the 
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other hand, if the wave-function at E is zero on the interface, the embedding 
potential blows up at the band edge like 



QmE oc 



1 



(28) 



E — Ei 



o 



Even though this one-dimensional model cannot be directly applied to the 
complicated bands of an fee transition metal, our eigenvalues behave like (j2*7|) 
or (|2*Hj) at band extrema, exhibiting the square root singularity on the em- 
bedding surface midway between atomic planes (figure 2b) where symmetry 
permits the band edge wavefunction to vanish. 

The interplay between the eigenvalues and the band-structure can be 
quite complicated. The tiny semi-circle in the eigenvalue figures, centred on 
—6.8 eV, varies in energy between the lowest Xi state and the maximum 
of the first Ai band. For each energy in this range there are two non-zero 
eigenvalues, corresponding to the two Ai Bloch states in the band structure. 
Similarly, there is a small eigenvalue semi-ellipse with energies between the 
minimum of the second A x band at E — —3.7 eV and its maximum at r 12 . It 
overlaps in energy with an outer eigenvalue starting and finishing at A = 0, 
in the range of energies for which there are again two Ai Bloch states. This 
topology shows that when there is more than one 0k,i at a particular energy 
with the same symmetry, there is not a one-to-one relationship with the ^k,? 
channel functions. In the cases we have just described, the two A x Bloch 
states are associated with one continuous band. This cannot have a one- 
to-one relationship with the two eigenvalues having the topology shown in 
figure 2. Rather, each Bloch state is a linear combination of all the channel 
functions at that energy, as given by (|24jl . Of course, in the case of a symme- 
try point like K = 0, the different bands at the same energy frequently have 
different symmetry, and then we can identify each band with one eigenvalue 
of E K - 

The ^-function peaks are associated with discrete energies in a symmetry 
gap, at E — —6.0 eV (Ai symmetry) and —0.5 eV (A5 symmetry). At these 
energies, the embedding potential has a singularity and its inverse S _1 has 
a zero eigenvalue, which from (J2J) means that there is a normal derivative 
satisfying 



In other words, there is a wave-function in the system, with an exponentially 
decaying envelope, with zero amplitude over the interface, and consequently 
infinite logarithmic derivative. This gives the 5-function peak in the eigen- 
value of Sk- This singularity only appears at an isolated energy in the limit 




(29) 
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of Sk being evaluated at a real energy with an infinitesimal positive imag- 
inary part, and as the associated eigenvector ipi(r s ) = there is no flux 
associated with it. 



4 Channel functions and scattering 

We now use our channel states to provide a simple proof of (JTJ, starting off 
by deriving a general scattering result containing SmS. Let us consider the 
geometry shown in figure 1 - we wish to find the full wave- function x(r) in 
region I, given an incoming incident wave V( r ) i n the substrate region II. 
This can be done in several ways starting from Green's theorem [HI EJ HQ , 
and it follows straightforwardly that x i n region I is given by 



G(r,r s )-p^(r s ) - J^-(r, r 5 )^(r 5 ) 
oris oris 



(30) 



where G is the Green function for the combined system of I and II; the surface 
normal in ()3())1 is taken to be outwards from region I. Now because G in fl3()(l 
is the outgoing Green function, we can use (J2J) to substitute for dG/dris as 
follows 

r)C r 

— (r,r 5 ) = -2 J s d 2 r' s G(r,r' s )Z(r' s ,r s ). (31) 
The incident wave-function, on the other hand, satisfies 

J^(rs) = -2 / s rf 2 r^*(r s ,r' s )^(r' 5 ), (32) 

the complex conjugation arising because this is an incoming wave. Substi- 
tuting (|3*T|) and (|3^j) into (|3TJj) gives the pleasing result 

X (r) = 2i [ d 2 r s [ d 2 r' s G{v, r s )3mS(r s , r' 5 )^(r' s ), r in I. (33) 
J s «* s 

At this stage we must be clear about the incident wave tp in (|3*3|) . This 
can be any wave in region II in the energy continuum (otherwise QmS in f!33|) 
vanishes), satisfying incoming boundary conditions. If region II is a semi- 
infinite bulk crystal, for example, ip is an arbitrary combination of the Bloch 
states approaching S, at fixed energy, together with exponentially decaying 
states. The incident wave is defined entirely by its value over S, ip(rs), which 
can take an arbitrary form, provided that it has some projection onto the 
flux-carrying states. 

We now apply (jSlIj) to a nanostructure connected to contacts over inter- 
faces Si and S r (figure 1). In dividing the system in this way, the nanos- 
tructure region also contains the atoms of the contacts which are perturbed 
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by the nanostructure. Region I consists of the whole structure to the right 
of Si, that is, the nanostructure plus the right-hand contact and lead. If 
ip(r) is the wave in the left-hand lead incident on the nanostructure, the full 
wave-function everywhere to the right is given by (jHHJ), with the double inte- 
gral evaluated over Si. We take the incident wave to correspond to the open 
channel function ip\ defined over the left-hand interface, normalised to unit 
flux over this interface. Then the wave-function over the right-hand interface 
is given by 

Xl (r r ) =2i\\ j d 2 nG{v r ,vi)4{vi), (34) 
Js l 

where A' is the non-zero channel eigenvalue (JHJ). We now expand x«( r r) 
in terms of the complete set of channel functions ipj{r r ) for the right-hand 
interface, 

Xi(rr) =^Uji>j(r r ), (35) 

3 

and normalising the right-hand open channel functions to unit flux we find, 
for open channels, 

% = 4iA^|/ d\ j dV r ^(r r )G(r r ,r0^i/). (36) 



Vt3 wi, 

" I Si JS; 



The transmission coefficient T i3 - is the transmitted flux in the jth channel 
corresponding to unit incident flux in the ith channel, and with our normal- 
isation this is given by 



Ta = \Uj\ 2 = 16(A^) 



d 2 n [ d 2 r r ^ r (r r )G(r r ,rO^(rz) 

S[ J S T 



2 



(37) 



This expression for the transmission coefficient is very simple, with the clear 
interpretation of G as a propagator from one channel to the other. 

Let us now sum T^j over all the open exit channels, to give the total 
transmitted flux for unit incident flux in channel i. Changing and simplifying 
the notation for spatial coordinates, this is given by 

T i3 = 16(A() 2 | dl r J dli J d3 r J d4,G(2,,l r ) 

x ( E (A5)V;(lrW(3r)| G*(3 r , 40^(40^(20- (38) 



open j 



K open j 



Now we have seen that SmS can be expanded in terms of the open channel 
eigenfunctions alone, and with unit flux normalisation (|T2|) becomes 



SmS(r s ,ry = -2 £ A 2 ^(r 5 )^(r' 5 ). (39) 



open i 
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Hence (|3*8*)l simplifies to 



E Tij = -8(A^) 2 [dl r Jd2 l Jd3 r /d4iG(2 z ,l r )9mE r (l r ,3 r ) 

open j 

xG*(3 r , 4,)Vi(4,)^(2,). (40) 
Now we sum over all incident channels, to give 

E = -8 J dl r J d2 t J d3 r J ^G^l^mE^lrA) 



open i, j 



X 



G*(3 r Ai)\ E (Ai) a ^(4,)V^(20 , (41) 



\ open ] 



and making use of (J39|) once again, we obtain the result 

E = 4 1 dl r / A | d3 r | d4,G(2,,l r )9taS r (lr,3 r ) 

open j, j 

xG*(3 r , A^QmE^i). (42) 

This gives the total transmission for unit flux in each of the incident channels. 

This result for total transmission is only useful if there is the same flux 
in each open incident channel. To understand this let us consider the case of 
semi-infinite interfaces, with unit flux in each of the incident Bloch functions 
at fixed K. Then the flux in the m'th channel function due to the flux in the 
i'th Bloch function is given in terms of the coefficients <2j m (j21j) by a* m ai m . 
So the total flux in the m'th channel due to all the open Bloch channels, 
each carrying unit flux, is given by J2i o.* m a im . But because a im is a unitary 
matrix, we have J2i a i m a im = 1, and unit flux in each Bloch function implies 
unit flux in each channel function. This is the result we require: with energy 
normalisation, the Bloch functions - hence the channel functions - all carry 
the same flux. We can then write the total transmission between left and 
right as 

T lr = 4 1 dl r J d2 t J d3 r J d4 / G(2i,l r .)SmE r (l n 3 r )G*(3 r) 4 / )SmEi(4i,2i) 
= 4Tr[Gj r 9mE r G*^mE,], (43) 

recovering (|T|). This result is manifestly independent of choice of channel 
functions. 

To summarise, we have shown how the eigenfunctions of the embedding 
potential, a generalised logarithmic derivative, may be used to define conduc- 
tion channels across a surface. These new channel functions are orthogonal, 
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which could instil distinct advantages in applications compared to Bloch 
states as conduction channels. Using the new channel functions we have pro- 
vided a simple derivation of a well-known result for the total transmission 
through a conductor system. 
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